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BACKGROUND: Mild brain hypothermia (31-34 °C) after 
neonatal hypoxia-ischemia (HI) improves neurodevelop- 
mental outcomes in human and animal neonates. Using an 
asphyxia model with neonatal mice treated with mild hypo- 
thermia after HI, we investigated whether 1 H nuclear mag- 
netic resonance (NMR) metabolomics of brain extracts could 
suggest biomarkers and distinguish different treatments and 
outcome groups. 

METHODS: At postnatal day 7 (P7), CD1 mice underwent right 
carotid artery occlusion, 30min of HI (8% oxygen), and 3.5 h 
of either hypothermia (31 °C) or normothermia (37 °C). Whole 
brains were frozen immediately after HI, immediately after 3.5 h 
of hypothermia or normothermia treatments, and 24 h later. 
Perchloric acid extractions of 36 metabolites were quantified by 
900 MHz ] H NMR spectroscopy. Multivariate analyses included 
principal component analyses (PCA) and a novel regression 
algorithm. Histological injury was quantified after HI at 5 d. 
RESULTS: PCA scores plots separated normothermia/HI ani- 
mals from hypothermia/HI and control animals, but more data 
are required for multivariate models to be predictive. Loadings 
plots identified 1 1 significant metabolites, whereas the regres- 
sion algorithm identified 6. Histological injury scores were sig- 
nificantly reduced by hypothermia. 

CONCLUSION: Different treatment and outcome groups are 
identifiable by ] H NMR metabolomics in a neonatal mouse 
model of mild hypothermia treatment of HI. 



Although mild therapeutic hypothermia after neonatal 
asphyxia improves outcomes and is now a standard of 
care, further insight into patient management is needed from 
studies on brain injury and recovery mechanisms, as well as 
into challenges in predicting outcomes. The ability to modify 
genes and their expression has made mouse animal models 
useful in elucidating pathological mechanisms. Publications 
already exist of murine investigations of hypothermia as a 
treatment for hypoxic-ischemic (HI) brain injury (1-4). 
Various murine-based investigations of other pathologies 



have found *H nuclear magnetic resonance (NMR) -based 
metabolomics to be useful (5-7). In a recent oxygen-glu- 
cose deprivation study with superfused respiring neonatal 
rat brain slices, we used l H NMR-based metabolomics of 
perchloric acid brain extracts to distinguish recovery groups 
belonging to different post-oxygen-glucose deprivation tem- 
perature treatments and to identify responsible metabolites 
(8). Large sample sizes, common in metabolomics studies of 
biofluids that can be obtained relatively easily, such as blood 
and urine, are often impractical and too costly in preclini- 
cal laboratory neuroscience studies. However, studies with 
smaller sample sizes can be useful if they produce experi- 
mental differences that are significantly large. In this study, 
l U NMR at 21.1 Tesla (900 MHz) was used for metabolomic 
analyses of whole brains, obtained from neonatal mice in 
an in vivo model of neonatal HI. If brain metabolomics can 
distinguish normothermic from hypothermic posttreatment 
data sets, future studies that add outcome data can address 
whether brain metabolomics can help guide therapy, predict 
outcomes, and become clinically useful. 

RESULTS 

Hypothermia Reduced Histological Injury 5 d After HI 

Figure 1 shows that mild hypothermia reduced brain injury 
in our model. The median injury score was 22.5 (range: 3-24) 
for the normothermia group as compared with 6 (range 3-20) 
for the mild hypothermia group. Therefore, mice that received 
hypothermia had lower injury scores as compared with mice 
treated with normothermia (P = 0.005). The overlap of the dis- 
tributions is small enough to consider using the terms "strong 
injury" and "mild injury," using 14 as the dividing score. Doing 
so leads to the conclusion that the hypothermia group had 
2.0 ± 1.3 mice (17 ± 11%) with strong injury, as compared with 
9.0 ±1.5 mice (75 ± 13%) with strong injury in the normother- 
mia group. A two -sample binomial proportion test to deter- 
mine if there is a statistical difference between 9/12 and 2/12 
resulted in a P value of 0.014. 



department of Anesthesia and Perioperative Care, University of California, San Francisco, San Francisco, California; department of Neurology, University of California, 
San Francisco, San Francisco, California; department of Epidemiology and Biostatistics, University of California, San Francisco, San Francisco, California; department of 
Pharmaceutical Chemistry, University of California, San Francisco, San Francisco, California; California Institute for Quantitative Biosciences (QB3), University of California, 
Berkeley, Berkeley, California; department of Pediatrics, University of California, San Francisco, San Francisco, California. Correspondence: Lawrence Litt (Larry.Litt@ucsf.edu) 

Received 1 7 November 201 2; accepted 1 April 201 3; advance online publication 26 June 201 3. doi:1 0.1 038/pr.201 3.88 



170 Pediatric RESEARCH 



Volume 74 | Number 2 | August 201 3 



Copyright © 201 3 International Pediatric Research Foundation, Inc. 



Neonate metabolomics in hypoxia-ischemia 



Articles 



Strong injury 




Mild injury 



Normothermia Hypothermia 

Figure 1 . Histological injury scores of mouse brains obtained after 5 d of 
recovery from hypoxia-ischemia (HI). Data at the left in red show scores 
for 1 2 mice kept normothermic after HI, whereas the data at the right in 
blue provide scores for 1 2 mice that were treated after HI with 3.5 h of mild 
hypothermia (31 °C). Median injury scores are shown by black horizontal 
bars. Mice that received hypothermia had lower injury scores as compared 
with mice treated with normothermia (P = 0.005). 

Nine Metabolites Were Identified in Pairwise f-Tests 
Figure 2 uses three vertical bar graphs to compare the 
quantifications of all 36 brain metabolites in five groups: 
the control group (control), the postnatal day 7 (P7) hypo- 
thermia group (P7 hypothermia), the P7 normothermia 
group (P7 normothermia), the P8 hypothermia group (P8 
hypothermia), and the P8 normothermia group (P8 normo- 
thermia). The detailed descriptions for each experimental 
group are provided in the Methods section. Figure 2a shows 
metabolites that had very significant metabolite changes, 
whereas Figure 2b shows metabolites to the left of the 4.67 
parts-per-million (ppm) water NMR peak, where metabo- 
lite signals are not detectable in vivo but are distinguish- 
able in this study because of the high magnetic field (21.1 
Tesla) and high spectral resolution (-0.002 ppm). Figure 2c 
shows the remaining metabolites. Eight of the nine metabo- 
lites identified in the t-tests were in Figure 2a: glutamine, 
acetate, succinate, aspartate, formate, N- acetyl -asp artyl- 
glutamate, isoleucine, and taurine. The ninth metabolite, 
fumarate, is in Figure 2b. Comparison results are summa- 
rized in Tables 1 and 2. 

Separate Clusters of Hypothermia and Normothermia Groups in 
Principal Component Analyses 

In Figure 3a, the scores plot for the plane of the first two 
principal component axes (R 2 = 0.55; Q 2 = 0.25) shows a 
clear separation of the HI group from all end of recovery 
groups, which overlap among themselves. In the corre- 
sponding loadings plot (Figure 3b), highlighting of the HI 
group's vector components shows that the HI group's sep- 
aration is due to increases in levels of alanine, ADP, cho- 
line, lactate, succinate, valine, Y- ammoDU tyric acid, and 



isoleucine, and decreases in levels of ATP, phosphocreatine 
(PCr), phosphocholine (PCho), malate, aspartate, taurine, 
and iV-acetyl-aspartate. 

Similarly, the scores plot in Figure 4a (R 2 = 0.563; Q 2 = 
0.002) for the ensemble of control and outcome groups shows 
a separation of the P8 normothermia group from the other 
groups, which are comingled. In the corresponding load- 
ings plot (Figure 4b), highlighting of the P8 normothermia 
group's vector components shows that the metabolite changes 
most responsible for the cluster separation were increases in 
glutamine, fumarate, succinate, isoleucine, iV-acetyl-aspartyl- 
glutamate, acetate, and formate, and decreases in taurine, his- 
tidine, malate, and ascorbate. 

Construction of the contribution plot (Figure 4c) began 
by calculating, one metabolite at a time, the average multi- 
variate space vector for each metabolite in the P8 normother- 
mia group. After bunching together all remaining treatment 
groups, the software similarly calculated each metabolite's 
average vector for the bunch. Green bars show, metabolite by 
metabolite, the magnitudes of vectors obtained by subtracting 
the average metabolite vectors in the P8 normothermia group 
from their partners in the bunched group. Numbers are plot- 
ted as positive or negative according to whether or not they 
are decreased or increased in the P8 normothermia group. 
Colored dots coordinate metabolite locations in Figure 4b 
with locations in Figure 4c. 

Regression Analyses: Good Linear Models Require Fewer Than 
Four Metabolites 

Results of the LI -penalized lasso regression analysis, per- 
formed with PCr as the outcome variable, are shown in 
Figure 5. The graph in Figure 5a shows the relationship 
between model predictive performance, as measured via 
mean squared error (MSE, j/-axis), and model complexity, as 
measured via log(X) (x-axis). The corresponding number of 
metabolites for a good linear model is given on the upper 
x-axis, at the top of the figure. Therefore, the PCr outcome 
results suggest that three to four variables are all that are 
needed to provide a best fit. 

Features of the graph that indicate very good success in the 
fitting include an initial drop in MSE, going from left to right, 
as the value of X starts to increase from ~0 (10~ 75 , because the 
log value is -7.5). The fits improve until the lowest point in 
the dip because unimportant metabolites are excluded and 
important metabolites increasingly show dominance. At a X 
value near 10~ 57 , the MSE has almost no change with fur- 
ther increases in A, suggesting that several statistical mod- 
els obtain approximately equal and minimal MSE. Finally, as 
X increases beyond 10~ 52 , one has ever-increasing values of 
MSE because not enough important variables are included. 
A companion plot (Figure 5b) shows coefficient trajecto- 
ries as successive metabolites enter the model, correspond- 
ing to relaxation of the LI penalty by decreasing X toward 
zero, going from left to right. Figure 5b shows that the top six 
l H NMR metabolites were ATP- ADP, PCho, lactate, glucose, 
valine, and glutamate. 
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Figure 2. Vertical bar graphs of all 1 H nuclear magnetic resonance (NMR) metabolites. Markings relate to results of pairwise comparison testing between 
groups: control vs. P7 hypothermia; control vs. P7 normothermia; control vs. P8 hypothermia; control vs. P8 normothermia; P7 hypothermia vs. P7 
normothermia; and P8 hypothermia vs. P8 normothermia. The five vertical bars with SEs for each metabolite correspond to protocol times when animal 
brains were taken. From left to right, the first bar refers to control animals (no hypoxia-ischemia (HI)), whereas the remaining four bars are for animals 
that underwent 30 min of HI followed immediately by 1 h of normoxia with the dam, followed by 3.5 h of normothermia or mild hypothermia, and then 
either no further recovery or 24 h of recovery. The vertical axis is for metabolite intensities normalized to dry powder weights and also to the 1 H signal 
intensity to the left of the DSS peak. Below each set of vertical bars is the multiplication factor used to raise each set's height to a more visible range, (a) 
Metabolites with statistically significant differences in between-treatment groups (*P< 0.05, + P< 0.01). (b) Metabolites with NMR peaks to the left of the 
water peak in Figure 1 c (*P < 0.05). (c) Metabolite quantifications not included in graphs a and b. DSS, 4,4-dimethyl-4-silapentane-1 -sulfonic acid; GABA, 
Y-aminobutyric acid; Gpcho, glycerophosphocholine; NAAG, A/-acetyl-aspartyl-glutamate; P, postnatal day; PCho, phosphocholine; PCr, phosphocreatine. 



Table 1 . Comparison results at 4.5 h after HI 



Metabolite 


Compared groups 


Tori 


Pvalue 


Aspartate 


P7 hypothermia vs. control 


I 


0.050 




P7 hypothermia vs. P7 normothermia 


i 


0.042 


Isoleucine 


Control vs. P7 normothermia 


T 


0.028 




P7 hypothermia vs. P7 normothermia 


T 


0.033 


Taurine 


P7 hypothermia vs. P7 normothermia 


T 


0.050 



HI, hypoxia-ischemia; R postnatal day. 



DISCUSSION 

The most important finding in this study is that although no 
metabolomics differences were apparent immediately after 
terminating the mild hypothermia treatment, at 24 h post- 
HI in the scores plot (Figure 4a), an impressive separation 
can be seen in the P8 normothermia group relative to the 
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hypothermia and control groups, which are indistinguishable 
from each other. 

The very clear separation of HI data sets from all others in the 
principal component analysis (PC A) scores plot (Figure 3a), 
which identified energy- related metabolites as primarily 
responsible, is more comforting than surprising, because HI 
brains were obtained at the time of maximal insult, before 
normoxia was restored, when conditions are physiologically 
very different from all other times. Other responsible variables 
included valine, a branched-chain amino acid that returns 
nitrogen to astrocytes in neuron-astrocyte glutamate recycling 
and that is also known to be increased during excitotoxicity 
(9). Astrocytic uptake of synaptic glutamate is followed by its 
amidation to glutamine. After glutamine is passed to presyn- 
aptic neurons for recycling into glutamate, the amino nitrogen 
taken from glutamine is returned to astrocytes via branched- 
chain amino acids (10). 
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In the HI and P8 normothermia groups, respectively, 
increased and decreased tricarboxylic acid cycle intermedi- 
ates succinate and malate can be linked to energy metabolism 
(11). The malate- aspartate shuttle is the primary mechanism 
for maintaining correct cytoplasmic and mitochondrial nico- 
tinamide adenine dinucleotide oxidized/nicotinamide adenine 
dinucleotide reduced ratios, whereas succinate contributes to 
complex II on the mitochondrial membrane. Lactate, alanine, 
and citrate can be transferred from astrocytes to neurons for 
metabolic support. The osmotic regulator taurine and the 
important neuron marker N- acetyl- aspartate are known to be 
decreased in neonatal HI. 

One report (12) using sensorimotor and water maze 
tests found persistent neurological deficits at 7 and 30 d 
after HI. Those results, together with our histology results, 
suggest that future HI murine studies with both outcome 
and metabolite data sets are likely to identify the metabo- 
lites most responsible for metabolomics differences and to 



Table 2. Comparison results at 24 h after HI 



Metabolite 


Compared groups 


Tori 


Pvalue 


Acetate 


P8 hypothermia vs. P8 normothermia 


I 


0.004 




Control vs. P8 normothermia 


i 


0.009 


Formate 


P8 hypothermia vs. P8 normothermia 


i 


0.014 




Control vs. P8 normothermia 


i 


0.028 


Fumarate 


Control vs. P8 normothermia 


i 


0.029 


Glutamine 


Control vs. P8 normothermia 


i 


0.014 


NAAG 


P8 hypothermia vs. P8 normothermia 


i 


0.045 


Taurine 


P8 hypothermia vs. P8 normothermia 


T 


0.050 


Succinate 


P8 hypothermia vs. P8 normothermia 


i 


0.049 



HI, hypoxia-ischemia; NAAG, /V-acetyl-aspartyl-glutamate; P, postnatal day. 



evaluate whether they are early predictors of long-term neu- 
rological outcomes. 

Both PCA scores plots (Figures 3a and 4a) had R 2 values 
>0.5, indicating that the two PCA statistical models explained 
>50% of each plot's total variance. However, for both PCA 
analyses, the Q 2 values were <0.25, which is less than the 
desired minimum value of 0.5. Therefore, although the find- 
ing of different clustering in multivariate space is significant, 
the model created from our data set was not sufficiently robust 
to be predictive. Clearly, we have too few data points for sus- 
taining a PCA cross-validation procedure that eliminates input 
data. In a clinical setting, one would need enough statistics to 
have predictive data sets, because these are used to alter man- 
agement, especially to try avoiding poor outcomes when they 
are predicted. 

The contribution plot (Figure 4c), which presents lengths 
of vectors in the multidimensional multivariate space, con- 
firms the clustering suggested in Figure 4b, which displays 
two-dimensional (2D) distances that are qualitative in a visual 
inspection of Figure 4b. Moreover, it confirms that there are 
no outliers and that identification of biomarker candidates is 
the same for both plots. 

The LI -penalized regression algorithm was very successful 
in identifying the metabolites and roles of the most important 
biomarkers. This method supplemented our PCA approach to 
biomarker identification, which is most successful in situations 
in which there are many multivariate data sets ("large AT'), and 
thus enough points on the scores plot for different validation 
tests of cluster separation. In PCA, the identity of metabolites 
responsible for separate clusters is found in the loadings plots. 

The LI -penalized regression identifies biomarkers in situa- 
tions for which it is not practical to have a very large N, which 
is the case in preclinical animal studies such as the current 
one and others, for which each point in a multivariate data 



2 

c\j 

O 0 
Q_ 

-2 
-4 




0 
PC 1 



0.2 



•Taurine 
• ATP 

PCrl I Pcho 
• Malate 
• ATP-ADP 



Glycine f a 

Phenylalanine* !j isti t din j; 

I Tryptophan 



I Glucose 



• Tyrosin 
I Glutathione »A 

I Myo-inositol 



IGpcho 



INTP 

I Aspartate 
I NAA 

•Threonine 
NAAG • |Cr 

• Acetate 
Glutamine | iGlutamate 
1 Formate 
• Fumarate 



0 
PC 1 



p| Alanine 



ADP# 

Choline! | Lactate 



Succinate 
GABA # %Va|jne 



I Isoleucine 
# Leucine 



Figure 3. Principal component analysis (PCA) scores and loadings plots for 18 data points, three from the end of HI, and 15 from the end of the recovery 
period, (a) Scores plot: horizontal and vertical axes are, respectively, for PCI and PC2. Colors: yellow, control; light green, P7 hypothermia; light blue, P8 
hypothermia; purple, P7 normothermia; reddish brown, P8 normothermia; black, HI. In a, black HI data are bunched to the right, separated from the 
outcome groups. R 2 values for PCI and PC2 are 0.386 and 0.1 67, respectively, (b) Loading plot: metabolites responsible for separations are shown in 
black bold print, next to brown dots; metabolites are shown in gray print next to green dots. GABA, y-aminobutyric acid; Gpcho, glycerophosphocholine; 
HI, hypoxia-ischemia; NAAG, A/-acetyl-aspartyl-glutamate; P, postnatal day; PCho, phosphocholine; PCr, phosphocreatine. 
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Figure 4. Principal component analysis (PCA) for the five outcome groups compared in Figure 2. Data were analyzed exactly as in Figure 3, after 
excluding the three HI data points, (a) Scores plot: data for the P8 normothermia group (red squares) are clustering in the right lower quadrant. 
R 2 values for PCI and PC2 are 0.243 and 0.1 75, respectively, (b) Loadings plot: metabolites contributing to the separation are highlighted in red and blue, 
(c) Contribution plot between the groups, as described in the text. Vertical green bars are distances in multivariate space between average metabolite 
vectors for each group. For both b and c, points in red and blue correspond to P8 hypothermia group metabolites that are increased and decreased, 
respectively. GABA, y-aminobutyric acid; Gpcho, glycerophosphocholine; HI, hypoxia-ischemia; NAAG, A/-acetyl-aspartyl-glutamate; P, postnatal day; 
PCho, phosphocholine; PCr, phosphocreatine. 



set costs thousands of dollars. For our small-N data set, the 
LI -penalized regression analysis produced easily recogniz- 
able minima in the MSE plots in Figure 5a,b. PCr is an out- 
come variable associated with a healthy state of high-energy 
phosphates, and the metabolites found by the algorithm are 
mechanistically relevant, with roles being apparent in Figure 
5b. All of the metabolites with negative coefficients occur in 
larger quantities when there is more injury and PCr deple- 
tion, whereas metabolites with positive coefficients occur 
in larger quantities when there is PCr preservation or aug- 
mentation. PCho increases can be produced by increased 
PCho synthesis via choline kinase, decreased activity of 
CTPxholine-phosphate cytidylyltransferase (EC 2.7.7.15), 
or breakdown of phosphatidylcholine (13). The significance 
of increased PCho being associated with greater energy 
preservation after HI in our neonatal mice is not totally 



clear, but one can speculate that the origin is the repair of 
damage or resuming of normal growth. The relative strength 
of metabolite roles is also apparent when reading Figure 5b 
from left to right. The ATP-ADP resonance peak is by far 
the most important, maintaining dominance as models with 
more metabolites are allowed. 

Although the NMR studies were performed in brain 
extracts, brains were removed after there was considerable 
time for significant metabolite exchanges between brain 
and blood. A recent search for biomarkers in blood taken 
<5 min after asphyxia in a nonhuman primate model of neo- 
natal asphyxia (14) used 2D gas chromatography with time- 
of-flight mass spectrometry and found significant changes 
in several blood metabolites, including succinate, lactate, 
glutamate, malate, and leucine. Although these metabolites 
do not freely traverse an intact blood-brain barrier, further 
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Figure 5. Plots produced by the LI -penalized lasso regression, (a) The vertical axis is for mean-squared errors (MSEs) from the regression's 
cross-validation procedure, plotted as a function of log(X), shown along lower horizontal axes. X Represents the tuning parameter. Numbers along the 
upper horizontal axis represent the average number of predictor variables, n, for linear models for different values of log(A). Going along a horizontal axis 
from right to left has A decreasing toward zero and n increasing toward 35. Vertical black lines bracket n values, where the model provides its best fits to 
the data. Red dots indicate average MSE values for all models resulting from the corresponding value of X, and vertical bars through the red dots show 
upper to lower MSE values, (b) Coefficient trajectories produced by the regression analysis (in a). The lower horizontal axis shows LI , the A-dependent 
sum of absolute values of (predicted - observed) differences. The upper horizontal axis indicates the number of predictor values shown in a. When going 
from left to right, A is decreasing toward zero. PCho, phosphocholine. 



studies of blood metabolites in our mouse model would be 
needed to assess their contributions to the brain metabolites 
that we quantified. 

What one would like to see, and which was not possible in 
this study, is metabolomics scores and loadings plots for later 
times after HI, ideally with one or more outcome scores associ- 
ated with each point in the multidimensional metabolite space. 
In an ideal situation, a metabolomics approach would have 
different clusters for different insult severities early after HI, 
different clusters for different treatment groups, and best out- 
come scores associated very strongly with one or more treat- 
ment groups. 

It is natural to ask when and how this study's complex 
NMR methods and statistical analyses might become clini- 
cally relevant and if research paths using in vivo magnetic 
resonance spectroscopy (MRS) are appropriate at this time. 
A major limit at this time is the spectral resolution of in vivo 
MRS at available magnetic fields. In 2003, the US Food and 
Drug Administration specified that risks to neonates younger 
than 1 mo have not been ruled out above 4 Tesla, nor have 
they been ruled out above 8 Tesla for older children and 
adults (15). Nevertheless, impressive technological advances 
in l W NMR metabolic profiling suggest a bright future. In a 
preclinical adult rodent study, 20 metabolites were accurately 
quantified in 14.1 Tesla in vivo l H MRS (16). In an in vivo 4.0 
Tesla human MRS clinical study of adults, the 2D J-resolved 
spectroscopy method was used to improve spectral resolu- 
tion (17). Moreover, a totally new MRS approach is evolving 
based on dynamic nuclear polarization spectroscopy, a tech- 
nology that is more complex and more expensive, but 35,000 



times more sensitive than existing systems at NMR detection 
of 13 C metabolites (18). 

In summary, this study demonstrated an important role for 
metabolomics in distinguishing different treatment and out- 
come groups, and it has identified potentially important bio- 
markers. We are encouraged that roles will ultimately emerge 
for multivariate metabolomics in helping clinical management. 

METHODS 

Animal Procedures: Induction of HI in Neonatal Mice and Mild 
Hypothermia 

All animal experimental procedures were approved by the University of 
California, San Francisco institutional animal care and use committee. 
For each HI experiment, a litter of 10 CD1 mice of either sex (Charles 
River Laboratories, Wilmington, MA) underwent the Vannucci proce- 
dure of neonatal HI (19-21) on P7. Briefly, each pup was placed on a 
37 °C temperature-controlled pad for general anesthesia from spontane- 
ous respirations of 3.0% isoflurane in oxygen. After the right common 
carotid artery was permanently ligated, pups recovered from anesthesia 
on a warming pad with littermates before being returned to the dam 
for 90min of feeding. HI was induced thereafter by placing the pups 
for 30min into chambers that were maintained at 37 °C and through 
which a mixture of 8% oxygen and 92% nitrogen flowed. During this 
period, one pups temperature was monitored continuously by a therm- 
istor temperature filament (outside diameter: 0.82 mm; probe type: 
RET-4; Physitemp Instruments, Clifton, NY) attached to the abdomi- 
nal skin. Temperatures of all other pups were checked regularly using 
a noncontact infrared laser thermometer. After HI, pups were placed 
with the dam for lh of recovery. They were then exposed to 3.5 h of 
either hypothermia (31 °C) or normothermia (37 °C). During hypo- 
thermia or normothermia, pup core temperatures were monitored as 
described above. Pups in the normothermia group were returned to 
the dam for 24 h after HI. Those in the hypothermia group underwent 
a slow rewarming (0.1 °C per minute) to 37 °C, the temperature of the 
normothermia group, and were then returned to the dam. 
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Figure 6. Using one nuclear magnetic resonance (NMR) tube with extracts from control mice, two different types of 1 H NMR spectra were obtained 
sequentially without moving the sample or changing static magnetic fields. A standard one-pulse } H NMR high-resolution 900 MHz spectrum 
(top tracing) is aligned above a ] H pJRES spectrum, the latter being a one-dimensional projection along the 1 H ppm axis of a two-dimensional J-resolved 
*H NMR plot, as discussed in the Methods. Notable differences in the pJRES spectra include its substantially flatter baseline and, for many metabolites, the 
consolidation of many resonance peak clusters into one peak. For the horizontal axis, parts-per-million (ppm) assignments relative to the DSS reference 
peak are from published chemical shifts. Spectral regions are (a) 0.8-2.9 ppm; (b) 2.9-4.2 ppm; (c) 5.6-8.8 ppm. Ala, alanine; Asp, aspartate; Cr, creatine; 
DSS, 4,4-dimethyl-4-silapentane-1 -sulfonic acid; GABA, Y-aminobutyric acid; Gin, glutamine; Glu, glutamate; Gly, glycine; GPCho, glycerophosphocholine; 
His, histidine; lie, isoleucine; Leu, leucine; NAA, A/-acetyl-aspartate; NAAG, A/-acetyl-aspartyl-glutamate; PCho, phosphocholine; PCr, phosphocreatine; Phe, 
phenylalanine; Thr, threonine; Trp, tryptophan; Tyr, tyrosine; Val, valine. 
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Figure 7. Expanded view of resonance changes for high-energy phosphate metabolites in six experimental groups, (a) ATP and ADP are shown for a 
spectral region generally free of peaks from nucleotide phosphates containing guanosine, cytosine, and thymidine, (b) The spectral region showing PCr 
and Cr being clearly separated. Cr, creatine; PCr, phosphocreatine. 



At each predetermined time point listed below, two pups, one at a 
time, were decapitated under isoflurane anesthesia, after which the 
heads were placed in ice. The two brains, after being removed rapidly 
and then snap-frozen in liquid nitrogen, were saved for subsequent 
perchloric acid extraction and NMR analysis of the mixed extracts. 
The time interval between decapitation and snap-freezing was always 
<60 s, which resulted in significant preservation of ATP and other 
high-energy phosphates that are known to disappear quickly in the 
brains of mice after decapitation (22,23). 

Using the notation P7 for postnatal day 7 pups and P8 for postna- 
tal day 8 pups, animals were assigned to the following six groups, from 
which brains were taken for NMR analysis: (i) The HI group: P7 brains 
obtained immediately at the end of HI, without reoxygenation; (ii) the 
control group (control): no HI, no hypothermia; (iii) the P7 hypother- 
mia group (P7 hypothermia): animals underwent the full HI protocol 
and hypothermia treatment, and brains were obtained immediately 
at the end of 3.5 h of hypothermia; (iv) The P7 normothermia group 
(P7 normothermia): animals underwent the full HI protocol and a 
3.5 h normothermia treatment, at the end of which their brains were 
obtained; (v) the P8 hypothermia group (P8 hypothermia): animals fol- 
lowed the same procedure as those in the P7 hypothermia group but 
underwent a slow rewarming before returning to the dam, with brains 
being obtained 24 h after HI; and (vi) the P8 normothermia group (P8 
normothermia): animals followed the same procedure as those in the 



P7 normothermia group but returned to the dam after treatment, with 
brains being obtained at 24 h after HI. For NMR analysis, a total of 36 
pups were used. Three NMR tubes were obtained from six mice for each 
of the six groups. 

Perchloric Acid Metabolite Extraction 

Each perchloric acid extraction was obtained from two whole brains 
of sibling P7 mice in the same experimental group. Perchloric acid 
extraction was done as in earlier studies (24-26). In brief, to obtain 
extracts for one NMR tube for a given time point, two frozen brains 
of sibling P7 mice in the same experimental group and time point 
were ground together under liquid nitrogen in a precooled mortar 
and then extracted using 30ml/g (dry mass) of ice-cold 12% perchlo- 
ric acid. After centrifugation, each supernatant was neutralized with 
KOH to pH 7.0-7.2 and then centrifuged again. The final supernatant 
was filtered through a 0.22 (im sterile MCE Fisherbrand syringe filter 
(Fisher Scientific, Pittsburgh, PA) and lyophilized for 24 h on VirTis 
BenchTop 2K lyophilizer (SP Industries, Gardiner, NY). 

NMR Data Acquisition and Analysis 

Before NMR analysis, lyophilized dry extracts were weighed and 
dissolved in 280 ul of D 2 0 containing 1 u.mol/1 4,4-dimethyl- 
4-silapentane-l -sulfonic acid, a reference compound for NMR 
determinations of chemical shifts and signal intensities. Samples 
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were neutralized to pD 7.0-7.2 with DC1 and NaOD before being 
loaded into 3 mm NMR tubes (Select Series S-3-900-7; Norell, 
Landisville, NJ). Spectra were acquired at the Central California 
900 MHz NMR Facility, which operates with National Institutes of 
Health support (GM68933) in the QB3 facilities at the University 
of California, Berkeley (Berkeley, CA). One-dimensional *H spectra 
were acquired on a Bruker Avance II 900 MHz (21.1 Tesla; Bruker, 
Billerica, MA) NMR spectrometer equipped with a 5 mm CPTCI 
cryoprobe (Bruker) using a 30° tip angle pulse with 128 K points 
over 4.65 s and with a recycle delay of 1.0 s, for a total interexperi- 
ment delay of 5.65 s. No solvent suppression was used. Spectra were 
processed without apodization or zero filling and were referenced to 
4,4-dimethyl-4-silapentane-l -sulfonic acid. The line widths of the 
methyl resonances of 4,4-dimethyl-4-silapentane-l -sulfonic acid 
typically ranged from 1.0 to 1.5 Hz (0.0011-0.0017 ppm). 

A check on *H resonance assignments for metabolites was done for 
each of the temperature conditions by performing 2D *H J-resolved 
NMR spectroscopy (2D JRES spectroscopy), using acquisition meth- 
ods described elsewhere (27). In a 2D JRES spectrum, units for the 
horizontal axis are ppm, and Hz units for the vertical axis are used 
to display distances of J-coupled peaks from each resonances center. 
Projection downward along the ppm axis consolidates the multiplets 
of each metabolite into a single peak, thereby greatly reducing the 
resonance peak overlap for different metabolites. The resulting pro- 
jection of the 2D JRES plot is known as a pJRES *H spectrum. 

Figure 6a-c show the full spectral regions, comparing one-pulse data 
(top tracings) with one-dimensional pJRES data (bottom tracings). 
Figure 7 shows expanded views of two small spectral regions contain- 
ing resonance peaks whose quantifications were important to the study. 
In Figure 7a, singlet resonance peaks from the adenosine moiety in 
ATP (8.49 ppm) and ADP (8.59 ppm) can be found. Quantification of 
the uncontaminated ADP was straightforward, but quantification of 
ATP involved deconvolution because of an overlapping component. 

*H NMR metabolites corresponding to different chemical shifts 
were identified and preprocessed with TopSpin 3.1 and AMIX soft- 
ware (Bruker) and quantified with iNMR (Nucleomatica, Molfetta, 
Italy). All chemical shift assignments for peaks and clusters were 
checked against chemical shift assignments published in careful stud- 
ies (28,29). Relaxation time corrections were not performed because 
the same NMR pulse sequence was used for all runs and the goal of 
the analysis was primarily to detect large differences in each metab- 
olite from its control value. Metabolite signal intensities were com- 
puted relative to 4,4-dimethyl-4-silapentane-l -sulfonic acid and then 
normalized both to the weight of the dry powder that was mixed with 
D 2 0 before NMR measurements and to the sum of quantifications of 
all metabolites shown in Figure 6. 

Brain Histology Assessment of Injury 

Although it is known that mild hypothermia is therapeutic in HI mouse 
models of neonatal asphyxia (1), data were obtained to validate and 
quantify this for the exact conditions of our experiment. Twenty- four 
additional P7 mice underwent the same HI protocol already described, 
with 12 mice assigned randomly to each of the normothermia and 
hypothermia groups. Animals were euthanized 5 d after HI, and brains 
were perfused via intracardiac injection of buffered (0.1mol/l phos- 
phate, pH 7.2) 10% formalin freshly prepared from paraformaldehyde. 
Coronal sections were cut through the forebrain at 50- urn intervals with 
a Vibratome (Ted Pella, Redding, CA). Alternate sections were stained 
with cresyl violet for morphology, or with Perl's stain to localize iron 
deposition. Brains were scored in a blinded manner as in earlier studies 
(20) for degree of injury, using a scale from 0 (no injury) to 24 (cystic 
infarct in multiple regions). Metabolomic data were not obtained. 

Statistical and Metabolomic Analyses 

Metabolite quantifications for all six experimental groups underwent 
within-group and between-group statistical analysis including pair- 
wise t- testing, PCA, and the LI -penalized lasso regression analyses. 

Pairwise comparison testing for differences in individual metabo- 
lites was performed via t- tests between groups treated with hypother- 
mia and normothermia, as well as between those groups and control 
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or HI groups. Therefore, there were 10 tests for each metabolite: P7 
hypothermia vs. P7 normothermia; P8 hypothermia vs. P8 normo- 
thermia; control vs. P7 hypothermia; control vs. P7 normothermia; 
control vs. P8 hypothermia; control vs. P8 normothermia; HI vs. P7 
hypothermia; HI vs. P7 normothermia; HI vs. P8 hypothermia; and 
HI vs. P8 normothermia. Multiplicity correction was effected using 
false- discovery rate methods (30). 

In the PCA, done using SIMCA P+ v. 13 software (Umetrics, San Jose, 
CA) as described previously (8,31), metabolite quantifications were cen- 
tered and scaled to SD units and analyzed as vectors in a multivariate 
space, with one point (or one data vector) in the space corresponding 
to a full set of metabolite quantifications. Two-dimensional scores and 
loadings plots were made using the first and second principal compo- 
nents (PCI and PC2), which correspond respectively to axes having 
the most and second-most variance for data vector projections. Scores 
plots were examined for separate groupings or clusterings of data from 
different treatment groups. Loadings plots, made from projecting indi- 
vidual metabolite vector components onto the plane of scores plots, 
were searched for the most important metabolites that contribute to 
clustering in the scores plots by examining metabolite components in 
the same quadrant as the clusters and maximally from the origin. Scores 
plots were assessed from two parameters, R 2 and Q 2 , these being known 
respectively as the "explained variance" and the "cross-validation vari- 
ance," or that fraction of the total variance that can be attributed to the 
vector components in the scores plot (32). 

Besides the scores and loadings plots, another PCA multivari- 
ate plot, the contribution plot, was used to compare metabolite dif- 
ferences between selected groups. This process starts with finding, 
metabolite by metabolite, each groups average multivariate space vec- 
tor for the chosen metabolite. For each metabolite, there is a vector 
subtraction of one groups metabolite vector from the others, followed 
by a calculation of the resulting vectors magnitude, which is then 
entered into the contribution plot. 

A novel regression algorithm, the LI penalized lasso (33,34), was 
employed as in our earlier studies to identify the most important metab- 
olites in linear models that explain changes in PCr, which we used as 
a surrogate for ATP. The lasso (least absolute selection and shrinkage 
operator) has recently emerged as a forefront regression technique for 
situations in which the number of covariates (here metabolites) is large 
as compared with the number of observations. As explained previously, 
the lassos tuning parameter, X, shrinks (reduces) the number of metabo- 
lites in the linear model as it increases upward from zero, initially elimi- 
nating noise variables, and ultimately eliminating the most important 
variables. The algorithm finds the optimal range for X (35). 

Median scores of brain injury for hypothermia and normothermia 
groups were calculated using Prism 5.0 (GraphPad Software, La Jolla, 
CA). The Mann- Whitney U- test was used to determine significance. 
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